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Ab inito calculations for the nuclear many-body problem make predictions for 
the density and isospin dependence of the nuclear equation-of-state (EOS) far 
away from the saturation point of nuclear matter. I compare predictions from 
microscopic and phenomenological approaches. Constraints on the EOS derived 
from heavy ion reactions, in particular from subthreshold kaon production, as 
well as constraints from neutron stars arc discussed. 



1. Introduction 

Heavy ion reactions provide the only possibility to reach nuclear matter 
densities beyond saturation density po ~ 0.16 frn -3 . Transport calculations 
indicate that in the low and intermediate energy range ~ 0.1 -j- 1 

AGeV nuclear densities between 2 -j- 3pn are accessible while the highest 
baryon densities (~ 8/30) wm probably be reached in the energy range 
of the future GSI facility FAIR between 20 ~ 30 AGeV. At even higher 
incident energies transparency sets in and the matter becomes less baryon 
rich due to the dominance of meson production. The isospin dependence of 
the nuclear forces which is at present only little constrained by data will 
be explored by the forthcoming radioactive beam facilities at FAIR/GSI , 
SPIRAL2/GANIL and RIA 2 . Since the knowledge of the nuclear equation- 
of-state (EOS) at supra-normal densities and extreme isospin is essential for 
our understanding of the nuclear forces as well as for astrophysical purposes, 
the determination of the EOS was already one of the primary goals when 
first relativistic heavy ion beams started to operate in the beginning of the 
80ties. In the following I briefly discuss the knowledge about the nuclear 

EOS at moderate densities and temperatures. For more details see e.g. Ref. 
3 



Models which make predictions on the nuclear EOS can roughly be di- 
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vided into three classes: phenomenological density Junctionals such as Gogny 
or Skyrme forces ^' ^' ^ and relativistic mean field (RMF) models 7 , effec- 
tive field theory (EFT) and ab initio approaches. In EFT a systematic ex- 
pansion of the EOS in powers of density, respectively the Fermi momentum 

& Q 

kp is performed. EFT can be based on density functional theory ■ or 
e.g. on chiral perturbation theory ^ Ab initio approaches are ased 
on high precision free space nucleon-nucleon interactions and the nuclear 
many-body problem is treated microscopically. Predictions for the nuclear 
EOS are essentially parameter free. Examples are variational calculations 
13 5 Brueckner-Hartree-Fock (BHF) ^ ^ or relativistic Dirac-Brueckner- 
Hartree-Fock (DBHF) ^ and Greens functions Monte-Carlo ap- 

proaches 20 j n follwoing I will mainly concentrate on the DBHF 
approach. 



2. The EOS from ab inito calculations 
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Fig. 1. Nuclear matter saturation points from relativistic (full symbols) and non- 
relativistic (open symbols) Brueckner-Hartree-Fock calculations based on different 
nucleon-nucleon forces. The diamonds show results from variational calculations. Shaded 
symbols denote calculations which include 3-body forces. The shaded area is the empir- 
ical region of saturation. Figure is taken from Ref. 



21 



In ab initio calculations based on many-body techniques one derives the 
energy functional from first principles, i.e. treating short-range and many- 
body correlations explicitcly. A typical example for a successful many-body 
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approach is Brucckncr theory . In the relativistic Brueckner approach the 
nucleon inside the medium is dressed by the self-energy S. The in- medium 
T-matrix which is obtained from the relativistic Bethe-Salpeter (BS) equa- 
tion plays the role of an effective two-body interaction which contains all 
short-range and many-body correlations of the ladder approximation. Solv- 
ing the BS-equation the Pauli principle is respected and intermediate scat- 
tering states are projected out of the Fermi sea. The summation of the 
T-matrix over the occupied states inside the Fermi sea yields finally the 
self-energy in Hartree-Fock approximation. This coupled set of equations 
states a self-consistency problem which has to be solved by iteration. 

In contrast to relativistic DBHF calculations which came up in the late 
80ties non-relativistic BHF theory has already almost half a century's his- 
tory. The first numerical calculations for nuclear matter were carried out by 

99 

Brucckncr and Gammcl in 1958 . Despite strong efforts invested in the 
development of improved solution techniques for the Bethe-Goldstone (BG) 
equation, the non-relativistic counterpart of the BS equation, it turned out 
that, although such calculations were able to describe the nuclear satura- 
tion mechanism qualitatively, they failed quantitatively. Systematic studies 
for a large number of NN interactions were always allocated on a so-called 
Coester-line in the E/A~p plane which does not meet the empirical region 
of saturation. In particular modern one-boson-exchange (OBE) potentials 
lead to strong over-binding and too large saturation densities where rela- 
tivistic calculations do a much better job. 

Fig. 1 compares the saturation points of nuclear matter obtained by 
relativistic Dirac-Brueckner-Hartree-Fock (DBHF) calculations using the 

no 

Bonn potentials ° as bare NN interactions to non-relativistic Brueckner- 
Hartree-Fock calculations for various NN interactions. The DBHF results 
are taken from Ref. ^ (BM) and more recent calculations based on im- 
proved techniques arc from Ref. ^' (Tubingen). Several reasons have been 
discussed in the literature in order to explain the success of the relativis- 
tic treatment (see e.g. discussion in Ref. ^). Three-body forces (3-BFs) 
have extensively been studied within non-relativistic BHF ^ and varia- 
tional calculations . Both approaches shwon in Fig. 2 are based on the 
latest AVi8 version of the Argonne potential. The variational results shown 
contain boost corrections (Sv) which account for relativistic kinematics and 
lead to additional repulsion . 

The contributions from 3-BFs are in total repulsive which makes the 
EOS harder and non-relativistic calculations come close to their relativis- 
tic counterparts. The same effect is observed in variational calculations io 
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Fig. 2. Predictions for the EOS of symmetric nuclear matter from microscopic ab initio 
calculations, i.e. relativistic DBHF ^ , non-relativistic BHF ^ and variational ^ cal- 
culations. For comparison also soft and hard Skyrme forces are shown. Figure is taken 
from Rcf. 21 . 



shown in Fig. 2. It is often argued that in non-relativistic treatments 3-BFs 
play in some sense an equivalent role as the dressing of the two-body interac- 
tion by in-medium spinors in Dirac phenomenology. Both mechanisms lead 
indeed to an effective density dependent two-body interaction V which is, 
however, of different origin. One class of 3-BFs involves virtual excitations 
of nucleon-antinucleon pairs. Such Z-graphs are in net repulsive and can be 
considered as a renormalization of the meson vertices and propagators. A 
second class of 3-BFs is related to the inclusion of explicit resonance degrees 
of freedom. The most important resonance is the A(1232) isobar which pro- 
vides at low and intermediate energies large part of the intermediate range 
attraction. 

Fig. 2 compares the equations of state from the different approaches: 

17 9*^ 1/1 

DBHF from Ref. ' based the Bonn A interaction* , BHF and varia- 
tional calculations The latter ones are based on the Argonne AVis po- 
tential and include 3-body forces. All the approaches use modern high pre- 
cision AW interactions and represent state of the art calculations. Two phe- 



a The high density behavior of the EOS obtained with different interaction, e.g. Bonn B 
or C is very similar. ^ 
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nomcnological Skyrme functionals which correspond to the limiting cases of 
a soft (K=200 MeV) and a hard (K=380 MeV) EOS are shown as well. In 
contrast to the Skyrme interaction where the high density behavior is fixed 
by the parameteres which determine the compression modulus, in micro- 
scopic approaches the compression modulus is only loosely connected to the 
curvature at saturation density. DBHF Bonn A has e.g. a compressibility 
of K=230 MeV. Below 3p both are not too far from the soft Skyrme EOS. 
The same is true for BHF including 3-body forces. 

When many-body calculations are performed, one has to keep in mind 
that elastic NN scattering data constrain the interaction only up to about 
400 MeV, which corresponds to the pion threshold. NN potentials differ 
essentially in the treatment of the short-range part. A model independent 
representation of the NN interaction can be obtained in EFT approaches 
where the unresolved short distance physics is replaced by simple contact 
terms. In the framework of chiral EFT the NN interaction has been com- 
puted up to N 3 LO 27 A n alternative approach which leads to sim- 
ilar results is based on renormalization group (RG) methods 28 . In the 
Mow k approach a low-momentum potential is derived from a given realis- 
tic NN potential by integrating out the high-momentum modes using RG 
methods. When applied to the nuclear many-body problem low momentum 
interactions do not require a full resummation of the Brueckner ladder dia- 
grams but can already be treated within second-order perturbation theory 

on 

. However, without repulsive three-body- forces isospin saturated nuclear 
matter was found to collapse. Including 3-BFs first promising results have 
been obtained with V\ ow k , however, nuclear saturation is not yet de- 
scribed quantitativley. 

2.1. EOS in symmetric and asymmetric nuclear matter 

Fig. 3 compares now the predictions for nuclear and neutron matter from 
microscopic many-body calculations - DBHF 10 and the 'best' variational 
calculation with 3-BFs and boost corrections - to phcnomcnological ap- 
proaches and to EFT. As typical examples for relativistic functionals we 
take NL3 ^ as one f the best RMF fits to the nuclear chart and a phe- 
nomenological density dependent RMF functional DD-TW from Rcf. 
ChPT+corr. is based on chiral pion-nucleon dynamics ^ including con- 
densate fields and fine tuning to finite nuclei. As expected the phenomeno- 
logical functionals agree well at and below saturation density where they 
are constrained by finite nuclei, but start to deviate substantially at supra- 
normal densities. In neutron matter the situation is even worse since the 
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Fig. 3. EOS in nuclear matter and neutron matter. BHF/DBHF and variational cal- 
culations are compared to phenomenological density functionals NL3 and DD-TW and 
ChPT+corr.. The left panel zooms the low density range. The Figure is taken from Ref. 
3 

isospin dependence of the phenomenological functionals is less constrained. 
The predictive power of such density functionals at supra-normal densities 
is restricted. Ab initio calculations predict throughout a soft EOS in the 
density range relevant for heavy ion reactions at intermediate and low en- 
ergies, i.e. up to about three times pq. There seems to be no way to obtain 
an EOS as stiff as the hard Skyrme force shown in Fig. 2 or NL3. Since the 
nn scattering lenght is large, neutron matter at subnuclear densities is less 
model dependent. The microscopic calculations (BHF/DBHF, variational) 
agree well and results are consistent with 'exact' Quantum-Monte-Carlo 
calculations . 

In isospin asymmetric matter the binding energy is a functional of the 
proton and neutron densities, characterized by the asymmetry parameter 
= Y n — Y p which is the difference of the neutron and proton fraction 
Yi — pi/ p ,i = n,p. The isospin dependence of the energy functional can 
be expanded in terms of (3 which leads to a parabolic dependence on the 
asymmetry parameter 



Fig. 4 compares the symmetry energy predicted from the DBHF and vari- 
ational calculations to that of the empirical density functionals already 
shown in Fig. 3 In addition the relativistic DD-prS RMF functional is in- 




(1) 
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eluded. Two Skyrme functionals, SkM* and the more recent Skyrme-Lyon 
force SkLya represent non-relativistic models. The left panel zooms the low 
density region while the right panel shows the high density behavior of 
-E sym . Remarkable is that most empirical models coincide around p ~ 0.6po 
where -E S ym — 24 MeV. This demonstrates that constraints from finite nu- 
clei are active for an average density slightly above half saturation density. 
However, the extrapolations to supra-normal densities diverge dramatically. 
This is crucial since the high density behavior of E sym is essential for the 
structure and the stability of neutron stars (see also the discussion in Sec. 
V.5). The microscopic models show a density dependence which can still be 
considered as asy-stiff. DBHF ° is thereby stiffcr than the variational rc- 
suits of Rcf. . The density dependence is generally more complex than in 
RMF theory, in particular at high densities where E syul shows a non-linear 
and more pronounced increase. Fig. 4 clearly demonstrates the necessity to 
constrain the symmetry energy at supra-normal densities with the help of 
heavy ion reactions. 

The hatched area in Fig. 4 displays the range of £ sym which has been ob- 
tained by constructing a density dependent RMF functional varying thereby 

no oo 

the linear asymmetry parameter 04 from 30 to 38 MeV . In Ref. 00 it 
was concluded that charge radii, in particular the skin thickness r n — r p 
in heavy nuclei constrains the allowed range of 04 to 32 -j- 36 MeV for 
relativistic functionals. 
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2.1.1. Effective nucleon masses 

The introduction of an effective mass is a common concept to character- 
ize the quasi-particle properties of a particle inside a strongly interacting 
medium. In nuclear physics exist different definitions of the effective nu- 
cleon mass which are often compared and sometimes even mixed up: the 
non-relativistic effective mass rn* NR and the relativistic Dirac mass m* D . 
These two definitions are based on different physical concepts. The non- 
relativistic mass parameterizes the momentum dependence of the single- 
particle potential. The relativistic Dirac mass is defined through the scalar 
part of the nucleon self-energy in the Dirac field equation which is absorbed 
into the effective mass m* D = M + Y^s(k, kp). The Dirac mass is a smooth 
function of the momentum. In contrast, the nonrelativistic effective mass 
- as a model independent result - shows a narrow enhancement near the 
Fermi surface due to an enhanced level density . For a recent review on 
this subject and experimental constraints on m* NR see Ref. ^ . 

While the Dirac mass is a genuine relativistic quantity the effective mass 
m* NR is determined by the single-particle energy 

(2) 

m* NR is a measure of the non- locality of the single-particle potential U (real 
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Fig. 5. The effective mass in isospin symmetric nuclear matter as a function of the 
momentum k at different densities determined from relativistic Brueckner calculations. 
Figure is taken from Ref. ^6 
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part) which can be due to non-localities in space, resulting in a momen- 
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turn dependence, or in time, resulting in an energy dependence. In order 
to clearly separate both effects, one has to distinguish further between the 
so-called k-mass and the E-mass . The spatial non- localities of U are 
mainly generated by exchange Fock terms and the resulting k-mass is a 
smooth function of the momentum. Non-localities in time are generated by 
Brueckner ladder correlations due to the scattering to intermediate states 
which are off-shell. These are mainly short-range correlations which gen- 
erate a strong momentum dependence with a characteristic enhancement 
of the E-mass slightly above the Fermi surface ' 1 . The effective 
mass defined by Eq. (2) contains both, non-localities in space and time and 
is given by the product of k-mass and E-mass .In Fig. 5 the nonrela- 
tivistic effective mass and the Dirac mass, both determined from DBHF 
calculations are shown as a function of momentum k at different Fermi 
momenta of kp — 1.07, 1.35, 1.7 fxn - . m* NR shows the typical peak struc- 
ture as a function of momentum around kp which is also seen in BHF 
calculations . The peak reflects the increase of the level density due to 
the vanishing imaginary part of the optical potential at kp which is also 
seen, e.g., in shell model calculations ■ . One has, however, to account 
for correlations beyond mean field or Hartree-Fock in order to reproduce 
this behavior. Fig. 6 compares the density dependence of the two effective 
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Fig. 6. Nonrelativistic and Dirac effective mass in isospin symmetric nuclear matter as 
a function of the density for various models. 



masses determined at kp. Both masses decrease with increasing density, the 
Dirac mass continously, while m^R starts to rise again at higher densities. 
Phenomenological density functionals (QHD-I, NL3, DD-TW) yield sys- 
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tcmatically smaller values of m* NR than the microscopic approaches. This 
reflects the lack of nonlocal contributions from short-range and many-body 
correlations in the mean field approaches. 



2.1.2. Proton-neutron mass splitting 

A heavily discussed topic is in the moment the proton-neutron mass split- 
ting in isospin asymmetric nuclear matter. This question is of importance 
for the forthcoming new generation of radioactive beam facilities which 
are devoted to the investigation of the isospin dependence of the nuclear 
forces at its extremes. However, presently the predictions for the isospin 
dependences differ substantially. BHF calculations ' 00 predict a proton- 
neutron mass splitting of m* NR n > m* NR . This stands in contrast to rela- 
tivistic mean-field (RMF) theory. When only a vector isovector p-meson is 
included Dirac phenomenology predicts equal masses m* D = m* Dp while 
the inclusion of the scalar isovector 5-meson, i.e. p+8, leads to m* D n < m* Dp 

on ' 

J . When the effective mass is derived from RMF theory, it shows the same 
behavior as the corresponding Dirac mass, namely m* NRn < m* NRv . 

on 

Conventional Skyrme forces, e.g. SkM*, lead to m* NR n < m* NRp oa while 
the more recent Skyrme-Lyon interactions (SkLya) predict the same mass 
splitting as RMF theory. The predictions from relativistic DBHF calcu- 
lations are in the literature still controversial. They depend strongly on 
approximation schemes and techniques used to determine the Lorentz and 
the isovector structure of the nucleon self-energy. Projection techniques are 
involved but more accurate and yield the same mass splitting as found 
in RMF theory when the 5 -meson is included, i.e. m* D n < m* Dp ^ 
Recently also the non-relativistic effective mass has been determined with 
the DBHF approach and here a reversed proton-neutron mass splitting was 
found, i.e. m* NR n > m* NR p Thus DBHF is in agreement with the results 
from nonrelativistic BHF calculations. 



2.1.3. Optical potentials 

The second important quantity related to the momentum dependence of the 
mean field is the optical nucleon-nucleus potential. At subnormal densities 
the optical potential U opt is constraint by proton-nucleus scattering data 
41 and at supra-normal densities constraints can be derived from heavy 
ion reactions, see Refs. ^ 43, 44 j n a re l a tivistic framework the optical 
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Schroedinger-equivalent nucleon potential (real part) is denned as 

One should thereby note that in the literature sometimes also an optical 
potential, given by the diff erence of the single-particle energies in medium 
and free space U = E- VM 2 + k 2 is used 42 which should be not mixed 
up with (3). In a relativistic framework momentum independent fields ^s,v 
(as e.g. in RMF theory) lead always to a linear energy dependence of U op t- 
As seen from Fig. 7 DBHF reproduces the empirical optical potential ^ ex- 
tracted from proton-nucleus scattering for nuclear matter at po reasonably 
well up to a laboratory energy of about 0.6-0.8 GeV. However, the saturat- 
ing behavior at large momenta cannot be reproduced by this calculations 
because of missing inelasticities, i.e. the excitation of isobar resonances 
above the pion threshold. When such continuum excitations are accounted 
for optical model caculations are able to describe nucleon-nucleus scattering 
data also at higher energies ^ . In heavy ion reactions at incident energies 
above 1 AGeV such a saturating behavior is required in order to reproduce 
transverse flow observables . One has then to rely on phenomenological 
approaches where the strength of the vector potential is artificially sup- 
pressed, e.g. by the introduction of additional form factors ^ or by energy 
dependent terms in the QHD Lagrangian ^ (D 3 C model in Fig. 7) . 

The isospin dependence, expressed by the isovector optical potential 
fiso = (^opt.n — [/ op t, P )/ (2/3) is much less constrained by data. The knowl- 
edge of this quantity is, however, of high importance for the forthcoming 
radioactive beam experiments. The right panel of Fig. 7 compares the pre- 
dictions from DBHF ^ and BHF ^ to the phenomenological Gogny and 
Skyrme (SkM* and SkLya) forces and a relativistic T — p approximation 
49 based on empirical NN scattering amplitudes At large momenta 
DBHF agrees with the tree- level results of Rcf. ^ . While the dependence 
of J7; so on the asymmetry parameter (3 is found to be rather weak ^> ^ ', the 
predicted energy and density dependences are quite different, in particular 
between the microscopic and the phenomenological approaches. The energy 
dependence of Ui so is very little constrained by data. The old analysis of 
optical potentials of scattering on charge asymmetric targets by Lane ^ is 
consistent with a decreasing potential as predicted by DBHF/BHF, while 
more recent analyses based on Dirac phenomenology ^ come to the oppo- 
site conclusions. RMF models show a linearly increasing energy dependence 
of Ui so (i-c quadratic in k) like SkLya, however generally with a smaller 
slope (see discussion in Ref. ° ). To clarify this question certainly more 
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Fig. 7. Nucleon optical potential in nuclear matter at pq. On the left side DBHF calcula- 
tions for symmetric nuclear matter from ^ an( j 17 are com p are( j to the phenomenological 

models NL3 and D 3 C ™ and to the p-A scattering analysis of "* . The right panel com- 

1 8 47 

pares the iso- vector optical potential from DBHF ° and BHF ' to phenomenological 
RMF ^ ( Gogny and Skyrme forces and to a relativistic T — p approximation ^ . 

experimental efforts are necessary. 

2.2. Probing the EOS by kaon production in heavy ion 
reactions 

With the start of the first relativistic heavy ion programs the hope was 
that particle production would provide a direct experimental access to the 
nuclear EOS 53 . It was expected that the compressional energy should be 
released into the creation of new particles, primarily pions, when the matter 
expands However, pions have large absorption cross sections and they 
turned out not to be suitable messengers of the compression phase. They 
undergo several absorption cycles through nucleon resonances and freeze 
out at final stages of the reaction and at low densities. Hence pions loose 
most of their knowledge on the compression phase and are not very sensitive 
probes for stiffness of the EOS. 

After pions turned out to fail as suitable messengers, K + mesons were 
suggested as promising tools to probe the nuclear EOS . At subthreshold 
energies K + mesons are produced in the high density phase and due to 
the absence of absorption reactions they have a long mean free path and 
leave the matter undistorted by strong final state interactions. Moreover, 
at subthreshold energies nucleons have to accumulate energy by multiple 
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scattering processes in order to overcome the threshold for kaon production 
and therefore these processes should be particularly sensitive to collective 
effects. Within the last decade the KaoS Collaboration has performed sys- 
tematic measurements of the K + production far below threshold, see Refs. 
55, 56, 57, 58, 59 Based on the new data situation, in Ref. ^ the question 
if valuable information on the nuclear EOS can be extracted has been re- 
visited and it has been shown that subthreshold K + production provides 
indeed a suitable and reliable tool for this purpose. These results have 
been confirmed by the Nantes group later on ^1. I n subsequent investiga- 
tions the stability of the EOS dependence has been proven, Refs. ^ . 
Excitation functions from KaoS ^ are shown in Fig. 8 and compared to 
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Fig. 8. Excitation function of the K"+ multiplicities in Au + Au and C + C reactions. 
RQMD 60 an d IQMD ^3 w ith in-medium kaon potential and using a hard/soft nuclear 
EOS are compared to data from the KaoS Collaboration ^ 



RQMD DU ' 11 and IQMD Dd calculations. As expected the EOS dependence 
is pronounced in the Au+Au system while the light C+C system serves as 
a calibration. The effects become even more evident when the ratio R of 
the kaon multiplicities obtained in Au+Au over C+C reactions (normalised 
to the corresponding mass numbers) is built Such a ratio has the 

advantage that possible uncertainties which might still exist in the theo- 
retical calculations should cancel out to large extent. This ratio is shown 
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in Fig. 9. Both, soft and hard EOS, show an increase of R with decreasing 
energy down to 1.0 AGeV. However, this increase is much less pronounced 
when the stiff EOS is employed. The comparison to the experimental data 
from KaoS where the increase of R is even more pronounced, strongly 
favours a soft equation of state. Fig. 9 demonstrates also the robustness 
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Fig. 9. Excitation function of the ratio R of K + multiplicities obtained in inclusive 
Au+Au over C+C reactions. RQMD ™ and IQMD ^3 calculations are compared to 
KaoS data ^ Figure is taken from . 



of this observable. Exploring the range of uncertainty in the corresponding 
transport calculations the stability of the conclusions drawn from this ob- 
servable has been demonstrated in Rcf. . This concerns elementary input, 
in particular the elementary production cross sections TV A; A A NYK + 
which are not constrained by data. 



2.3. Constraints from neutron stars 

Measurements of "extreme" values, like large masses or radii, huge lumi- 
nosities etc. as provided by compact stars offer good opportunities to gain 
deeper insight into the physics of matter under extreme conditions. There 
has been substantial progress in recent time from the astrophysical side. 

The most spectacular observation was probably the recent measure- 
ment on PSR J0751+1807, a millisecond pulsar in a binary system 
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with a helium white dwarf secondary, which implies a pulsar mass of 
2.1 ± 0.2 Mq with la (2a) confidence. Therefore, a reliable EOS has 

to describe neutron star (NS) masses of at least 1.9 M© (la) in a strong, 
or 1.6 M Q (2a) in a weak interpretation. This condition limits the softness 
of EOS in NS matter. One might therefore be worried about an appar- 
ent contradiction between the constraints derived from neutron stars and 
those from heavy ion reactions. While heavy ion reactions favor a soft EOS, 
PSR J0751+1807 requires a stiff EOS. The corresponding constraints are, 
however, complementary rather than contradictory. Intermediate energy 
heavy-ion reactions, e.g. subthreshold kaon production, constrains the EOS 
at densities up to 2 -f- 3 po while the maximum NS mass is more sensitive 
to the high density behaviour of the EOS. Combining the two constraints 
implies that the EOS should be soft at moderate densities and stiff at high 
densities. Such a behaviour is predicted by microscopic many-body calcula- 
tions (see Fig. 2). DBHF, BHF or variational calculations, typically, lead to 
maximum NS masses between 2.1 — 2.3 Mq and are therefore in accordance 
with PSR J0751+1807, see Ref. 65 . 

There exist several other constraints on the nuclear EOS which can be 
derived from observations of compact stars, see e.g. Refs. 65 - 66 - 67 . Among 
these, the most promising one is the Direct Urea (DU) process which is 
essentially driven by the proton fraction inside the NS DU processes, 
e.g. the neutron /3-decay n — > p + e~ + D e , are very efficient regarding their 
neutrino production, even in superfluid NM and cool NSs too fast to 

be in accordance with data from thermally observable NSs. Therefore, one 
can suppose that no DU processes should occur below the upper mass limit 

for "typical" NSs, i.e. Mrju > 1-5 Mq (1.35 Mq in a weak interpretation). 

71 

These limits come from a population synthesis of young, nearby NSs 
and masses of NS binaries ^4 

3. Summary 

The status of theoretical models which make predictions for the EOS can 
roughly be summarized as follows: phenomenological density functionals 
such as Skyrme, Gogny or relativistic mean field models provide high preci- 
sion fits to the nuclear chart but extrapolations to supra-normal densities or 
the limits of stability are highly uncertain. A more controlled way provide 
effective field theory approaches which became quite popular in recent time. 
Effective chiral field theory allows e.g. a systematic generation of two- and 
many-body nuclear forces. However, these approaches are low momentum 
expansions and when applied to the nuclear many-body problem, low den- 
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sity expansions. Ab initio calculations for the nuclear many-body problem 
such as variational or Brueckner calculations have reached a high degree 
of sophistication and can serve as guidelines for the extrapolation to the 
regimes of high density and/or large isospin asymmetry. Possible future de- 
vellopments are to base such calculations on modern EFT potentials and 
to achieve a more consistent treatment of two- and three-body forces. 

If one intends to constrain these models by nuclear reactions one has 
to account for the reaction dynamics by semi-classical transport models of 
a Boltzmann or molecular dynamics type. Suitable observables which have 
been found to be sensitive on the nuclear EOS are directed and elliptic 
collective flow pattern and particle production, in particular kaon produc- 
tion, at higher energies. Heavy ion data suggest that the EOS of symmetric 
nuclear matter shows a soft behavior in the density regime between one 
to about three times nuclear saturation density, which is consistent with 
the predictions from many-body calculations. Conclusions on the EOS are, 
however, complicated by the interplay between the density and the mo- 
mentum dependence of the nuclear mean field. Data which constrain the 
isospin dependence of the mean field are still scare. Promising observables 
are isospin diffusion, iso-scaling of intermediate mass fragments and particle 
ratios (tt + /tt~ and eventually K + / K° '^). Here the situation will certainly 
improve when the forthcoming radioactive beam facilities will be operating. 
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